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We investigate the cosmic evolution of the linear bias in the framework of a flat FLRW spacetime. 
We consider metric perturbations in the Newtonian gauge, including Hubble scale effects. Making 
the following assumptions, (i) scale independent current epoch bias feo, (ii) equal accelerations be- 
tween tracers and matter, (iii) unimportant halo merging effects (which is quite accurate for z < 3), 
we analytically derive the scale dependent bias evolution. The identified scale dependence is only 
due to Hubble scale evolution GR effects, while other scale dependence contributions are ignored. 
We find that up to galaxy cluster scales the fluctuations of the metric do not introduce a significant 
scale dependence in the linear bias. Our bias evolution model is then used to derive a connection 
between the matter growth index 7 and the observable value of the tracer power spectrum nor- 
malization 0-8(2;). We show how this connection can be used as an observational test of General 
Relativity on extragalactic scales. 
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1. INTRODUCTION 

The distribution of matter on large scales, based on dif- 
ferent extragalactic objects, can provide important con- 
straints on models of cosmic structure formation. How- 
ever, a serious problem that hampers such a straight- 
forward approach is our limited knowledge of how lumi- 
nous matter traces the underlying mass distribution. In 
particular, the concept of the so called biasing between 
different classes of extragalactic objects and the back- 
round matter distribution was put forward by Kaiser 
i| and Bardeen et al. Q in order to explain the higher 
amplitude of the 2-point correlation function of clusters 
of galaxies with respect to that of galaxies themselves. 

Such a biasing is statistical in nature; with galaxies 
and clusters being identified as high peaks of an under- 
lying, initially Gaussian, random density field, while it is 
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assumed to be linear and scale- independent"'^. Formally, 
the linear bias factor, 6, is defined as the ratio of the ex- 
tragalactic mass tracer fluctuations, (5tr, to those of the 
underlying mass, Sm'- 



b5„ 



(1.1) 



Since the two-point correlation function in a continuous 
density field is defined as ^(r) = {6{x)6{x + r)), one can 
write the bias factor as the square root of the ratio of 
the two-point correlation function of the tracers to the 
underlying mass: 



6r 



1/2 



(1.2) 



in which case one considers the large-scale correlation 
function, ie., scales corresponding roughly to the so- 
called halo-halo term of the dark matter (hereafter DM) 
halo correlation function (see for example 4]). Further- 
more, since the variance of a density field, smoothed 



^ However, in the non-linear scales of clustering, ie., mostly below 
5 Mpc, a scale dependence of the bias has been observed 

(see for example Q). 
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at some scale R, is the correlation function at zero lag 
(<T^ = ^_r(0) — (^^(x))) one can also write the linear 
bias factor as the ratio of the variances of the tracer and 
underlying mass density fields, smoothed at some linear 
scale, traditionally taken to be 8 h^^ Mpc (at which scale 
the variance is of order unity): 

6=^. (1.3) 

We refer the reader to a thorough discussion of the meth- 
ods to estimate the tracer correlation function and also 
of the corresponding determination of the tracer bias val- 
ues, at different redshifts, that appears in Papageorgiou 
et al. (Q; and references therein). 

The bias factor may have many dependencies; even 
assuming that it is scale independent, it necessarily de- 
pends on the type of the mass tracer as well as on the 
epoch z, since the fluctuations evolve with time as gravity 
draws together galaxies and mass. It is evident, therefore, 
that the bias factor should also depend on the different 
dark energy models (hereafter DE), including those of 
modified gravity Q- It is the redshift evolution of bias, 
b{z), which is very important in order to relate obser- 
vations with models of structure formation and it has 
been shown to be a monotonically increasing function of 
redshift 

In the literature there are two basic families of ana- 
lytic bias evolution models. The first, called the galaxy 
merging bias model, is based on the Press-Schechter Q 
formalism, on the peak-background split J3] and on the 
spherical collapse model @, and reproduces relatively 
well the results of numerical simulations, although dif- 
ferences have been found especially at the high and low 
DM halo mass range. These differences have lead to mod- 
ifications of the original model to include the effects of 
ellipsoidal collapse [Ifll ; to new values of the bias model 
parameters [ll|; to new forms of the bias model fitting 
function [l2| or even to a non-Markovian extension of the 
excursion set theory [T3| . 

The second family of bias evolution models assumes 
a continuous mass-tracer fiuctuation field, proportional 
to that of the underlying mass, and the tracers act as 
"test particles". In this context, the hydrodynamic equa- 
tions of motion and linear perturbation theory are used. 
An original suggestion, named galaxy conserving bias 
model used the continuity equation and the assumption 
that tracers and underlying mass share the same veloc- 
ity field ^M^, while the bias evolution is provided 
by the solution of a 1st order differential equation as: 
b(z) = 1 -I- (6o — with bo the bias factor at the 

present time and D{z) the growing mode of density per- 
turbations. However, this bias model suffers from two 
fundamental problems: the unbiased problem ie., the fact 
that an unbiased set of tracers at the current epoch re- 
mains always unbiased in the past, and the low redshift 
problem ie., the fact that this model represents correctly 
the bias evolution only at relatively low redshifts z<^ 0.5 
[i8i]. Note that [l^ has extended this model to also in- 



clude an evolving mass tracer population in a ACDM 
cosmology. 

An attempt to derive a bias evolution model free of the 
above mentioned problems, utilized all three hydrody- 
namical equations of motion, linear perturbation theory 
and the fact that mass-tracers and underlying mass share 
the same gravity field, but not necessarily the same ve- 
locity field and that the linear bias is scale independent. 
This resulted in a second order differential equation in 
6, the solution of which provided the evolution of the 
linear and scale- independent bias (see @, and [2]|). 
We would like to stress here that the provided solutions 
apply to cosmological models, within the framework of 
general relativity. 

In the context of a scale-independent bias factor, an 
extension of the previous model, valid for all DE and 
modified gravity cosmologies, was recently proposed by 
Basilakos, Plionis & Pouri [7|. This extension provides 
a tool, using the evolution of bias, to put constraints 
on those DE models which adhere to general relativity, 
as well as to investigate whether the DE reflects on the 
nature of gravity ( "geometrical dark energy" ) . 

Overall, the scope of the present article is (a) to ex- 
tend the original Basilakos et al. [3| bias solution, by 
taking into account possible contributions from the met- 
ric fluctuations, and (b) to propose new tools that can 
be used in order to test the validity of general relativity 
on cosmological scales. 

The structure of our paper is as follows. The basic the- 
oretical elements of the problem are presented in section 
II, where we introduce [for a spatially flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry] the ba- 
sic cosmological equations. The issue related with the 
linear bias is discussed in section III. In this section we 
also present the general bias solution in the framework, 
by taking into account metric perturbations for the New- 
tonian gauge. Finally, the main conclusions are summa- 
rized in section V. 



2. SCALE DEPENDENT MATTER AND 
TRACER DENSITY PERTURBATIONS 

Let us derive the basic equations that govern the evo- 
lution of the mass density contrast as well as of the ex- 
tragalactic tracers, modeled here as a presurreless fluid 
(Pm = Ptr = 0). Note that the perturbed FLRW space- 
time in the Newtonian gauge is given by 

ds^ = -(1 + 2$)rfi2 + (1 - 2^)a^{t)dx^ (2.1) 

where $ is the Newtonian potential, a{t) is the scale fac- 
tor (normalized to unity at the present time) and dx^ is 
the flat spatial metric. In the current paper we assume a 
slowly varying gravitational potential $. Note that con- 
sidering the unperturbed spacetime one can easily derive 
the background equations 
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SttG 



{Pm + Pde) , 



p + 3H{p+pde) = 



(2.2) 



(2.3) 



In the above set of differential equations, an over-dot de- 
notes derivative with respect to time, p„i and pde, are 
the matter and dark energy densities with p = p„i + pde, 
H = a/a is the Hubble parameter, whereas p^e = WdePde, 
corresponds to the pressure assuming non-clustering dark 
energy. Note that for Wde{z) = — 1 we recover the con- 
cordance ACDM model. 



A. Matter density perturbations 

In this section, we discuss the basic equation which 
governs the evolution of the matter perturbations up 
to horizon scales and within the framework of any DE 
model. Following the notations of Dent et al. [22[ the 
perturbed (anisotropic stress-free) equations in the New- 
tonian gauge take the form 



$ = -4:H<^> + 8TTGpdeWde<^ 



5m = 3$ + ^w/,, 



Vf.rn = 

with constraint equations 

k 



3i/(i/$ + $) + — $ = -47rGp™<5„ 



(2.4) 

(2.5) 
(2.6) 

(2.7) 



(i/$ + $) = -4TrGpmVf,m (2.8) 

where f/,m = —VmO. (vm is the velocity potential for mat- 
ter). In this context, the combination of the relativistic 
equations (|2.4p - (|2.8p obtains the basic differential equa- 
tion 



fc2 

2HSrn + - . 



(2.9) 



A solution of the above equation provides the evolution 
of the matter fluctuations in the linear regime. 

On the other hand, the linear matter overdensity Sm = 
Spm/Pm is written as a function of the gravitational po- 
tential $ and the background variables as follows [23] : 



inGpmSn 



(2.10) 



which is the usual Poisson equation. In this context, in- 
serting the Poisson equation p. lip into eq. (|2.9|) we derive 
the well known scale independent equation 



2HS,n - 4:TTGp,nS„i = 



(2.12) 



Now, for any type of dark energy an efficient parametriza- 
tion of the matter perturbations is based on the growth 
rate of clustering [2J| 



/o(«) 



dlna 



(2.13) 



where 7 is the so called growth index (see Refs. [25 
ill) and n„(a) = ilma-^ / E"^ (a) ."^ Integr ating eq.^J^ 
we obtain an approximate solution of eq. (|2.12p which is 
valid for any type of dark energy'^: 



Sm{a) = aexp 



° dx 

- (f^7„(^)-i) 



(2.14) 



where is the scale factor of the universe at which 
the matter component dominates the cosmic fluid (here 
we use Ui ~ 10~^). Following standard lines we have 
Sm{a) oc D{a), where D{a) is the linear growing mode, 
usually scaled to unity at the present epoch D{a) = 
Smio.) / 5m{^) ■ It is interesting to mention that measur- 
ing the growth index could provide an efficient way to 
discriminate between modified gravity models and DE 
models which adhere to general relativity. Indeed it was 
theoretically shown that for DE models inside general 
relativity the growth index 7 is well fit by 7gr « 6/11 
(see [H,!!!). 

For the benefit of the reader we point here that for 
the traditional A cosmology it has been found, by some 
of us [22, 35], that the linear matter fiuctuation field 
starts to become scale-dependent, due to metric pertur- 
bations in Newtonian gauge, on scales larger than about 
- 50-100/i-iMpc (fc < 0.01-0.02/iMpc-i). Therefore, 
for large scales we have to use the generalized Poisson 
equation (|2.10p which is valid up to horizon scales. No- 
tice that on dimensional grounds we may approximate 
the quantity 3i?$ in Eq. as 3H<^ ~ 3H^<^ (see 

also [22,]). This is justified on the basis of Eq. 12.41 since 
[given also that AnGpde = 0{H^)] the only timescale that 
determines the evolution of $ is the Hubble scale H. It is 
therefore a good approximation to assume that $ ~ 
Using the latter condition and inserting Eq. (|2.10p into 
Ea. (l2.9p one can easily find that 



2HSm - 47rG, 



effPm^m 



Sm = 



(2.15) 



where G denoting Newton's gravitational constant. In 
the sub-Hubble (small scale) approximation (fs- ^ H^) 
equation (|2.10l) takes the form 



AnGprnSr, 



(2.11) 



^ Qm is the density parameter at the present time and E{a) = 
H{a)/HQ is the normalized Hubble function. For the usual A 
cosmology we have E{a) = {flrnO,~^ -I- 1 — ^^m)^^'^- 

^ Since the pure matter universe (Einstein de-Sitter) has the solu- 
tion of 5es m = o.: normalize our DE models to get 5 
at large redshifts, which should hold due to the dominance of the 
non-relativistic matter component. 
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where 



and 



Geffia,k) 



G 



l + efe(a,fc) 



(2.16) 



(2.17) 



For many DE models, it is convenient to study the 
growth evohition in terms of the expansion scale a rather 
than t. If we change the variables from i to a (^ = 
aH-^) then the time evolution of the mass density con- 
trast (see Eq. (|2.15|) ) takes the following form 



A{a\ 



d6m 
da 



B{a,k)6 
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where 



and 



B{a,k) = 



Aia) 



d\nE 
da 



3 
a 



2a^E^{a) 



[{l + ^k{a,k)] 



(2.18) 



(2.19) 



(2.20) 



We would like to end this section with a discussion 
on the evolution of the scale dependent growth rate of 
clustering. Obviously, it becomes important to construct 
a scale-dependent parametrization that is analogous to 
Eq. p.lSp and solves (approximately) Eq. ()2.18|) for aU 
scales k. In order to construct such a parametrization 
we focus on the matter dominated era when most of the 
growth occurs and express Cfc(a, k) as 



^k{a,k) = 



ac^k^ 



(2.21) 



In this context. Dent et al. 22] proposed that the scale- 
dependent growth rate /(a, k) may be expressed in terms 
of the scale- independent growth rate /o(a) in the form: 



f{a,k) = 



dlnSmia, k) 



Ma) 



dlna 1 + ^fe(a, k) 

Thus from Eq. (|2.22p we simply get 

nz,{x) - 1 



(5m (a, k) = aexp 



" dx 

X 



l+a(a:,fc) 



l + efc(a,fc) ■ 
(2.22) 



(2.23) 



Due to Sm{a,k) cx D(a,k) the normalized growth factor 
becomes 



D{a,k) 



Srnja, k) _ 5m{z,k) 
S„,{l,k) ~ (5„(0,fc) 



(2.24) 



where a = (1 -1- z) ^. 

B. Tracer density perturbations 

Now we use the same formalism as before but for the 
tracers. We would like to spell out what are the basic 



assumptions here (see also [36|). First of all we con- 
sider that the mass tracer population is conserved with 
time i.e. that the effects of hydrodynamics (merging, 
feedback mechanisms etc) do not significantly alter the 
population mean. However, the effects of merging have 
been phenomenologically modeled and it has been found 
(using N-body simulations) that they are important only 
for 2.5 — 3 [2l| as far as the bias factor is concerned. 
In what follows we also treat tracers with ptr = which 
implies that galaxies (or clusters of galaxies) are colli- 
sionless. Thus the tracer density evolves as 



ptr + 3ffptr = 



(2.25) 



Due to the same gravity field the corresponding equations 
(|2.4p and (|2.7I) are also valid here. On the other hand we 
have 



Vf.tr 

(i/$ + $) 



3$+ „^ 

-AnGptrVf^tr 



(2.26) 
(2.27) 
(2.28) 



where w/.tr = —VtrO- and ftr is the velocity potential of 
the tracers (in general different with that of matter t;tr 7^ 
I'm), Ptr is the tracer density, pm is the mass density and 
<& is the gravitational potential. 

Since the tracers and the underlying matter share 
the same gravitational field, this implies that the gen- 
eralized Poisson equation (I2.10p remains practically the 
same. In other words we could have different velocity 
fields [vtr ^ Vm) but the corresponding accelerations 
(wtr = v) are the same. Again, by taking the approxi- 
mation 3i7$ ~ 3i/2$ and using Eqs. ([223), (EH), (Ell) 
and Eqs. (|2.26|) - (l2.28p . we can obtain after some algebra 
the time evolution equation for the tracer fluctuation field 







or 



da^ 



Aia) 



da 



B{a,k)Srn=0 



(2.29) 



(2.30) 



C. Gauge Dependence 

The Newtonian gauge used in the above calculations 
is physically interesting because it corresponds to a time 
slicing of isotropic expansion. However, the matter den- 
sity perturbation (5,„(t, k) is a gauge dependent quantity 
and therefore it is important to clarify how our results 
change in alternative gauges, and what their connection 
is with observable gauge invariant quantities. 

Scalar metric perturbations around a spatially fiat 
background can be written in the following general form 
341 
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where a and r are the conformal cosmic expansion scale 
factor and the conformal cosmic time; "|" denotes the 
background three-dimensional covariant derivative. The 
corresponding perturbed energy-momentum tensor Tj^ 
has the form 

T^a = /3m(l + 5m) , 
-^0 ^ PmU\i , 

T''j = -prn^\ij , (2.32) 

where pm is the unperturbed pressureless matter density; 
U and E determine velocity perturbation and anisotropic 
shear perturbation. 

Gauge choices simplify the above expressions by setting 
various quantities to 0. For example we have = $ = 
B = in the Synchronous Gauge (SG), B = E = 
in the Newtonian Gauge (NG) and U ^ B = in the 
Conioving Time-orthogonal Gauge (GTG). 

In the special case of the synchronous gauge, which 
corresponds to a time slicing obtained by the matter local 
rest frame everywhere in space (the free falling observer 
frame), the line element of the perturbed spacetime is 
given by 

ds^ = a^{T)[-dT^ + {5ij + hij)dx'dx^] (2.33) 

It is straightforward to derive the growth equation for 
6m = 5pm/pm in & matter dominated universe in the 
synchronous gauge to obtain [H, [3^ 

5f^ + 2H5^„^ - ^TrGpmS?,? = (2.34) 

This growth equation is exact in the synchronous gauge 
in the case of matter domination and involves no scale 
dependence as in the case of equation ()2.15|) of the Newto- 
nian gauge. This scale independence is an artifact of the 
particular time slicing of the synchronous gauge which is 
a good approximation on small scales but is unable to 
capture the horizon scale effects modifying the growth 
function on large scales. 

Nevertheless equations (|2.15[) and (|2.34p clearly agree 
on small scales where ^fe ^ 0. Therefore, for larger scales 
{k < O.OlhA'Ipc) the question that arises is the following: 
What is the proper gauge to use when comparing with 
observations? 

This question has been addressed in Ref. [s^l where 
a gauge invariant observable replacement was obtained 
for dm- This observable (S°f''(t, fc) involves the matter 
density perturbation Sm(t,k) corrected for redshift dis- 
tortions due to peculiar velocities and gravitational po- 
tential. It also includes volume and position corrections. 
The final expression however is complicated and makes 



(2.31) 

I 

the theoretical predictions based on it not e asy to im- 
plement and manipulate. However, in Ref. [3l| it was 
pointed out that the Newtonian gauge matter perturba- 
tion is a good approximation to the observable gauge 
invariant perturbation (5^*(i, k) even on very large scales 
(comparable to the Hubble scale). This result will also 
be justified in the remaining part of this section. 

The cosmological perturbations evolution is well de- 
scribed by the gauge-invariant (GI) approach, pioneered 
by Bardeen[3^. This approach may be used to identify 
physical observables as gauge invariant quantities (e.g., 
Refs.jH, [13). A gauge-invariant matter density pertur- 
bation may be constructed as [H, Issj 

Sg'^ = 6ra + 3-iU ~B) , (2.35) 
a 

S'^/^ coincides with the density perturbation 5^ in the 
Synchronous Gauge for the pressureless matter system. 
Thus 6^^^ corresponds to the density perturbation rela- 
tive to the observers everywhere comoving with the mat- 
ter. These free falling observers do not experience the 
isotropic expanding background of the universe because 
the peculiar velocity of matter is distinct from the Hub- 
ble fiow. Thus has physical significance only for 
perturbations on scales small compared to the Hubble 
scale. 

In addition to S^'^ , it is straightforward to construct 
alternative gauge invariant quantities related to the mat- 
ter overdensity as evaluated in different gauges. Such 
a gauge-invariant variable closely related to the matter 
overdensity in the Newtonian gauge is of the form [32| - 

Sg'^ = Sm + -{B-E)= S„, - 3-{B ~ E) . (2.36) 
p a 

and has important advantages over Sg^^ discussed in 
what follows. S^^^ coincides with the density pertur- 
bation (5^*^ in the Newtonian Gauge (NG), in which 
B = E = 0. 

It is also straightforward to construct two gauge- 
invariant scalar potentials (f) and ip, which reduce to the 
gravitational potential in the Newtonian limit: [33 |: 

cj>^^-^{B^E) , 
a 

^ = ^, + -^[{B-E)a] . (2.37) 
a dr 

The gauge invariant gravitational potential (f> obeys the 
Poisson eauation[3^ IsSj sourced by S^^^ with no appear- 
ance of the Hubble scale: 

= -fc20 = AnGpa^dg^'^. (2.38) 
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where k is the (comoving) wavenumber of the Fourier 
mode. The Poisson equation is vahd only for scales small 
compared to the Hubble radius 1/H while on scales larger 
than the Hubble scale the growth of matter density per- 
turbations is frozen. Hence, 6^^^ can not be regarded 
as the observable matter density perturbation on scales 
comparable to the Hubble scale. Therefore, the observ- 
able density perturbation on both the small-scale and the 
large-scale modes can not be described directly by S^^^ 
even though it is a gauge-invariant quantity. 

The other gauge invariant perturbation S^^^ has some 
important attractive features with respect to observabil- 
ity, not shared by 6^^^ . These are summarized as follows: 

• It reduces to the Newtonian gauge perturbation 
6^'~^ , i.e. it corresponds to a frame which respects 
the isotropic expansion of the universe and is there- 
fore more appropriate for description of large scale 
perturbations. This reduction also simplifies the 
calculation of this perturbation. 

• It drives a scale dependent modification of the Pois- 
son equation for the gauge invariant potential 0. 
Indeed, the time-tir ne p art of the linearized Ein- 
stein equation gives [U [s^l 



V 



2^ 3- + ' 
o a 



2 cGIN 



(2.39) 



Thus, the anticipated scale dependence on scales 
comparable to the Hubble scale is picked up by the 
perturbation S^^^ . 

• It is gauge invariant as anticipated for any observ- 
able quantity. 

Thus, the gauge invariant S^^^ and the Newtonian 
gauge variable 6^'~^ to which it reduces, constitute an 
attractive choice for making theoretical calculations to 
obtain the gravitational potential and the matter den- 
sity perturbation that can be directly compared with 
observations on large scales. However, these theoreti- 
cally obtained quantities need to also be corrected for 
bias, redshift distortions (due to gravitational potential 
and peculiar velocities), lensing magnification and vol- 
ume distortion |30i. 



3. THE EVOLUTION OF BIAS 

Clearly, due to the metric perturbations, equations 
(I2.15P and (|2.29p involve a scale k dependence of bias in 
contrast to the small scale approximate equation (j2.12p 
which is scale- invariant. Therefore, due to Eq. ()l.l|) one 
would expect that the bias factor must inherit a similar 
dependence to that of the density fluctuations namely 
b = h{z,k). Within this framework, we can distinguish 
three possible bias evolution cases: 

Case 1: Tracers and Mass .share same velocity field: 
Here we use the assumption of Tegmark fc Peebles (Tcj 



(see also [3), that the tracers and the underlying mass 
distribution share the same velocity field. Using the 
latter and Eqs. (f231) . (|2?26| we have 







(3.1) 



Now since we assume linear biasing, i.e. Eq. ()l.l|) . we 
obtain: 



where y — b — 1 and 5m oc D. An integration of Eq. 
provides: 



b{z,k) = l + y{z,k) = 1 + 



bo -I 
D{z,k) 



(3.3) 



where bo is the bias factor at the present time. 

This model is known to suffer from the so-called unbi- 
ased and the low redshift problems, by which the bias of 
mass tracers which either obey bo < 1 or are located at 
relative large redshifts, z > 0.5, cannot be modeled by 
Eq.dSSl) 

Case 2: Tracers and Mass share same acceleration 
field: Now we consider that both the tracers and the un- 
derlying mass distribution share the same gravitational 
field but different velocity fields [S]. Inserting Eg. dl.ip 
into Ea. (|2.15p and using simultaneously Eg. (12. 291) . we 
obtain a second order differential equation which de- 
scribes the evolution of the linear bias factor, b, between 
the background matter and the mass-tracer fluctuation 
field: 

ySm + 2((5™ + HSm)y + AirGcsPmSmy = , (3.4) 

where y = b — 1. ^ Transforming equation p.4p from 
t to a, we simply derive the evolution equation of the 
function y{a,k) [where y{a,k) = b{a,k) — 1] which has 
some similarity with the form of ea. (|2.18p as expected. 
Indeed, we have: 



fy 



Aia) 



2/(a,fc) 



dy 
da 



+ B{a,k)y = 0. (3.5) 



In Basilakos & Plionis (y, |20[ , we have provided an ap- 
proximate solution of Ea. p.4p . using f{z, k) = fo{z) ^ 1 
(which is valid at relatively large redshifts) , for cosmolog- 
ical models in the framework of general relativity, which 
contain quintessence (or phantom) dark energy. Here 
our aim is to provide an exact bias solution taking into 
account metric perturbations, namely b(z, fc), for all pos- 
sible dark energy cosmologies (for the case with no metric 
perturbations see 0). 



^ The current theoretical approach does not treat the possibiUty of 
having interactions in the dark sector. Also discussions beyond 
the linear biasing regime can be found in [37| (and references 
therein). 
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Inserting now y{a, k) = g{a)/D{a, k) into Ea. p.5|) and 
using simultaneously equation (j2.18p and the first equal- 
ity of equation ()2.22|) . we obtain: 



, A, -.dg 



a general solution of which is 



9W 



C1+C2 



da 



a^E{a) 



(3.6) 



(3.7) 



where Ci and C2 are the integration constants. Utilizing 
now a = (1 + z)~\ & = 2/ + 1 = [g/D) + 1 and Eq. ((X71) . 
we finally obtain the functional form which provides our 
general solution for all possible types of DE models, as: 



where 



b{z,k) = 1 + 



J{z) 



1 



C-2 



J(z) 



D{z,k) D{z,k) 



(1 + x)dx 



E{x) 



(3.8) 



(3.9) 



An extension of the above model to include the effects 
of halo merging processes, which introduces one fur- 
ther component in Ea. p.8|) . has been phenomenologically 
modeled in and it was found, using cosmological N- 

body simulations, that such effects are important only for 
2.5 — 3. Therefore, in the light of currently available 
" growth of structure" data (which reach z ^ 1; WiggleZ 
[38j), the merging term in the bias evolution model has 
been neglected. 

Notice that the dependence of our bias evolution model 
on the different cosmologies enters through the different 
behavior of D{z, k), which is affected by 7 (see equations 
[1231 [2;21, and oiE{z) = H{z)/Ho. Since different halo 
masses result in different values of Bq, one should expect 
that the constants of integration Ci = 60 — 1 and C2 
should be functions of the mass of dark matter halos (see 
[2]|). assuming that the extragalactic mass tracers are 
hosted by a DM halo of a given mass. It is interesting 
to mention here that our bias model, similarly to most 
others proposed in the literature, relate a mass tracer, be- 
ing a galaxy, an AGN or a cluster of galaxies, with a host 
dark matter halo within which the mass tracer forms and 
evolves. The models themselves follow the linear evolu- 
tion of the host halo and not the internal evolution of the 
astrophysical processes of the tracer. Thus the assump- 
tion is that the effects of nonlinear gravity and hydro- 
dynamics (merging, feedback mechanisms, etc.) can be 
ignored in the linear-regime (see [ill3). 

Comparing our solution of Eq. (|3.8p with that of Case 
1, b{z, k) — 1 + {bo — 1)/D{z, k), it becomes evident that 
the latter misses one of the two components of the full 
solution simply because the assumption of equal veloci- 
ties leads to a first order homogeneous differential equa- 
tion for the bias (13.21) . While the assumption of equal 
accelerations (with utr Vm) involves the full linear per- 
turbation analysis and thus it produces a second order 
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FIG. 1: The scale-dependent bias z-evolution {upper panel) 
at galaxy cluster scales, k ~ 0.05/iMpc~^ (open symbols), 
while the scale-independent prediction is shown as the solid 
line. In the lower panel we present the fractional difference 
with respect to the scale-independent bias model. Note that 
we use flm ~ 0.273, 'Wde{z) — —1, cr8,m{0) = 0.81 and 7 = 
0.55. 



homogeneous differential equation (j3.5p . independent so- 
lutions of which are 1/D{z,k) and J{z)/D{z,k). 

The further component in the bias solution, provided 
by the above model, solves the the known unbiased and 
the low redshift problems, by which the Tegmark & Pee- 
bles [1^ {Case 1) model suffers. 

Case 3: Tracers and Mass do not share same velocity 
field: Here we obtain a similar to Case 2 bias model but 
using the nomenclature of Tegrmark and Peebles [l^. 
We now drop the main assumption used in Case 1, that 
the mass-tracers and the underlying mass distribution 
share the same velocity fields, by allowing some sort 
of relation between the two (matter and mass-tracers) 
velocity fields. We obtain again the corresponding equa- 
tion p.2p . starting only from the continuity equations 
()2.5l I2.26P and introducing an additional time-dependent 
term, w/.tr(a) — w^.m(a), which we associate with the 
effects of the different velocity fields. We also use the 
same notation, as in our original formulation, that the 
tracers and the underlying mass distribution share the 
same gravity (accelerations) field. Then we obtain: 



Str - S,n = b/,tr(a) - W/,m(a)] 



(3.10) 



^^ = ^b/,t.(«)-«/,™(«)] (3.11) 



dt a 
a general solution of which is 



y{a,k)5m{a,k) = Ci + 



k'^da 
a^H{a) 



[vf,tr{a) - w/,m(a)] 



Thus the evolution of bias becomes: 

5o-l /(z,fc) 



b{z,k) = 1 + 



D{z,k) D{z,k) 



(3.12) 
(3.13) 
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where 



I{z,k) 



{1 + x) 



,0 E{.) ^^(oTfc)^!^^'-^^^-"^-^^)]'^^ 

(3.14) 

and hQ — \ = Ci/i5m(0, k). Obviously, if w/^tr = ^'/.m then 
Eq.dSHSl) boils down to Eg. ip^)) as it should. Although 
we do not have a fundamental theory to model the time- 
dependent 'y/,tr(a) — w/,m(a) function, it appears to de- 
pend on the Hubble constant as well as on the scale 
due to the fact that the quantity I{z, k) of Ea. (l3.14p has 
to be unit-less. With the above in mind, we thus observe 
that from Eqs-dS^, dSl]) and ([XTi)) we obtain: 



(3.15) 



C2-ffo'5m(a = l,fc) 

VfMa) - VfMa) = p 



implying the following scaling of the velocity potentials 
with the scale factor: 



vtria) - Vmia) 



l,k) 



ak'^ 



(3.16) 



In general the above difference between the velocity po- 
tentials could have had a different form with respect to 
the scale factor, but it must always be oc H^/k^ in order 
for the I{z, k) integral to be unit-less. 

Finally, in Figure 1 (upper panel) we compare the 
scale-dependent (symbols) and scale-independent (solid 
curve) bias evolution b{z, k) at galaxy cluster scales, k ~ 
0.05/iMpc-i (r c± 20/i-iMpc with M ~ 5 x IO^/i-IMq). 
As it is expected, the biasing is a monotonically increas- 
ing function of redshift. In the lower panel we show the 
fractional difference between the two bias evolution re- 
sults and it becomes quite evident that the fluctuations 
of the FLRW metric do not affect the bias evolution. 



4. TESTING GRAVITY AT COSMOLOGICAL 
SCALES 

Nowadays, the issue of testing the validity of general 
relativity at cosmological scales is considered one of the 
most fundamental and challenging problems on the inter- 
face uniting Astronomy, Cosmology and Particle Physics 
and indeed in the last decade there have been theoret- 
ical debates among cosmologists regarding the methods 
that researchers have to develop in order to make this 
achievement. Briefly, it is interesting to mention that 
measuring the growth index could provide an efficient 
way to discriminate between modified gravity models 
and DE models which adhere to general relativity. In- 
deed it was theoretically shown that for DE models in- 
side general relativity the growth index 7 is well fitted by 
7gr ~ 6/11 (see i28i],j2S]). Notice, that in the case of the 
braneworld model of Dvali, Gabadadze & Porrati [111 we 
have 7 ~ 11/16 (see also [23|), while for the f{R) gr avity 
models we have 7(2) « 0.41 - 0.21z for Vim = 0.273 So||- 

Recently, it has been proposed (see for example |41| ) 
that an efficient avenue to constrain the 7 parameter is by 




M^=5xl0i3h-iMQ ^ 

I I I I I I 



FIG. 2: The predicted growth history of the Universe for 
different flat cosmological models (upper panel: Mh = 6 x 
10^^/i"U'/q, middle panel: Mh = 2.7 x 10^^/i"^Mq) and bot- 
tom panel: Mh = 5 x 10^^ Mq) and their fractional dif- 
ference with respect to the ACDM model (see insert panels). 
The cosmological models shown are: f(R) (dashed hne) with 
7(2) = 0.41 - 0.212, concordance ACDM (solid hne) with 
7 = 0.55 and DGP (dot-dashed) with 7 = 0.68. Note that we 
use ftm = 0.273, as,m(0) = 0.81. 



determining observationally the redshift-dependent lin- 
ear growth of perturbations. Other methods have also 
been proposed in the literature, such as redshift space 
distortions in the galaxy power spectrum, weak lensing 
and the growth rate of massive galaxy clusters (see for 
example 42,| and references therein) . Indeed, it has been 
shown [431 that a good test to discriminate among "clas- 
sical" DE models and modified gravity models is to com- 
pare the expected combination parameter of the growth 
rate of structure, f(z), and the redshift-dependent rms 
fluctuations of the linear density field, crg(z), with that 
measured observationally (from large redshift surveys, 
like the WiggleZ; [H, [l^ and references therein). 

Armed, with our bias evolution model it is straightfor- 
ward to obtain theoretically a model-independent way of 
expressing the parameter combination f(z)(7g_tr(z)- Since 
the metric fluctuations do not signiflcantly affect the 
evolution of bias and thus the formation of large scale 
structures at the scales of interest, we utilize: £,k(z, k) = 
0. Therefore, using the solution of the bias evolution 
ea. (|1.3p and eq. (|2.22p we find that the growth history of 
the Universe is given by: 

f(z)a8Mz)iMh,z) - ni(z)b(Mh,z)asMz) (4.1) 

where as,m(z) = as,m{0)D(z) and b(Mh,z) is given by 
ea (|3.8l) . Our proposed gravity test consists in comparing 
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the above expectation with observationally determined 
values of /(z)(T8,tr(-2)) estimated by using existing or fu- 
ture redshift catalogs of extragalactic mass tracers (large 
red galaxies, optical or X-ray QSO, clusters of galaxies, 
etc). Note, that such data are already available in the lit- 
erature for the case of bright emission- line galaxies [sHi • 
Evidently, the essential cosmological parameters that 
enter in the theoretical expectation of Ea. (l4.ip are: 
^nnWde{z),as,m{0) and 7. Note however that: 

• there is only a weak dependence of 7 on w{z), as has 
been found in Linder & Cahn [28|, which implies 
that one can separate the background expansion 
history, E{z) = H{z)/Hq, constrained by a large 
body of cosmological data (SNIa, BAO, CMB), 
from the fluctuation growth history, given by 7, 
and 

• the value of cr8,T?i(0) remains relatively constant 
(<78,m(0) S [0.77,0.81]) for a range of dark energy 
equations of state (ACDM, quintessence, CPL), as 
shown by a recent analysis of SDSS Large Red 
Galaxies (H. 

In particular, the aim of our proposed method is to con- 
strain, for a given expansion history, the value of 7, and 
test whether there are deviations from the GR expecta- 
tions. In order to visualize the redshift and 7 dependence 
of the /cTs.tr, we compare in Figure 2 for different masses 
of dark matter halos, three flat cosmological models, in 
which we impose — 0.273 and crg^m (0) = 0.81. In 
particular, we consider the following cases: 

(a) the /(i?) model with 7(2) = 0.41 - 0.21z (dashed 
line) , 

(b) the concordance ACDM (7 = 0.55, solid line), and 

(c) the DGP model with 7 = 0.68 (dot-dashed hne). 
The inset panels of Figure 2 show the relative differ- 
ence of the f{R) or DGP model fas parameter combina- 
tion with respect to that of the ACDM. Interestingly, the 
f{R) models show the largest deviations (> 10%) at the 
lower redshift end (z < 0.5), while the DGP model shows 
large deviations also at the higher redshift end. It should 
be mentioned that qualitatively and quantitatively the 
relative differences among the models are quite similar, 
independently of the DM halo mass used, a fact which 
indicates that any extragalactic mass tracer can be used 
with the same efficiency to perform this cosmological test. 

We will elaborate on the details of our proposed 
method in a forthcoming paper, but it is important to 
realize that the parameters (C2, &o) of our bias model de- 
pend: (a) on the characteristic DM halo mass, within 
which the mass tracer is embedded (see Appendix A), 
and (b) on the values of fi™ and cr8.m(0) (see for exam- 
ple, Q). In any case we expect that the variation of 
the bias model (60, C2) parameters, within a physically 
acceptable range of Qm,o's,m(fi) values, should be quite 
small. For example, the fact that 6n[cx l/aH.m.(0)] with 
o-8,m(0) e [0.77,0.81] (as indicated in [13) implies that 
60 remains mostly unaffected as far as its dependence on 
o'8.m(0) is concerned. 



5. CONCLUSIONS 

In the current work we provide a general bias evolu- 
tion model, based on linear perturbation theory, which 
takes into account also metric fluctuations. We find that 
the metric fluctuations do not affect the evolution of bias 
and thus the formation of large scale structures. We ar- 
gue that the evolution of the mass-tracer fluctuations, 
quantified by their variance on some smoothed scale, can 
be used to test the validity of general relativity on cosmo- 
logical scales. We would like to remind the reader that in 
Basilakos et al. paper we have derived the evolution 
of bias within the context of scale-independent bias. The 
combination of the latter and current works provides a 
complete investigation of the linear bias evolution issue 
in cosmological studies. We show in our current work 
that the use of the combination parameter of the growth 
rate of structure and the rms fluctuations of the linear 
density field could provide an efficient avenue to discrimi- 
nate among "geometrical" (modified gravity) dark energy 
models and those that adhere to general relativity. 

It is however important to spell out clearly which are 
the basic assumptions of our model, which are common 
also to many other bias models in the literature: (a) Hub- 
ble scale GR effects taken into account in the fluctuation 
growth (b) the Newtonian gauge approach is employed 
(c) the mass tracers and the underlying mass share the 
same gravity field but different velocity fields, in contrast 
to the bias model proposed by Tegmark & Peebles [l^ 
in which they proposed the they share the same velocity 
field, (d) the biasing is linear on the scales of interest, and 
(e) that each DM halo is populated by one extragalactic 
mass tracer, which is an assumption that enters, at the 
present development of our model, only in the compari- 
son of our model with observational bias data and not in 
the derivation of its functional form. Finally, we assume 
unimportant halo merging effects which is quite accurate 
for z < 3. 
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Appendix A: Parametrizing the Bias Evolution 
Model using N-body Simulations 



Our analytical solution Eq. p.Sp gives a family of dark 
matter (DM) halo bias curves with two unknown parame- 
ters {bo, C2), which depend on the halo mass as well as on 
the cosmological background (see [lO] and the Appendix 
in 0). One can determine the behavior of the linear bias 
factor as a function redshift and halo mass, evaluating 
these constants using, for example, N-body simulations. 
To this end we use the results of the high resolution, 
collisionless, WMPA7 ACDM simulation of [1]. Here we 
only present the basic information regarding this N-body 
simulation. 

The simulation is a random realization of the concor- 
dance ACDM model (4^J with a volume of a 500 Mpc 
cube and 512^ particles. The adopted cosmological pa- 
rameters are: fim = 0.273, ft a = 1 — flm, h — 0.704 and 
as = 0.81, while the particle mass is 7.07 x 10^°h~^MQ 
comparable to the mass of a single galaxy. The sim- 
ulation was performed with an updated version of the 
parallel Tree-SPH code GADGET2 46]. 

Furthermore, the details of the method used to iden- 
tify the DM halos and estimate their bias, as a function 
of redshift, with respect to the underlying mass distribu- 
tion have been presented elsewhere (eg., [45j and refer- 
ences therein) . We only mention here that the DM halos 
are defined using a Friends of Friends algorithm with a 
linking length I = 0.17{n)~^^^ , where {n) is the mean 
particle density and that the DM halo bias is estimated 



by measuring the ratio of the variance of the DM halo 
fluctuations to that of the underlying DM in spheres of 
8 Mpc radius, while its uncertainty is based on the 
bootstrap re-sampling technique. 

For our present analysis we use 10 different redshift 
snapshots of the N-body simulation, spanning the range: 
< z < 5, while the DM halo catalogs are determined 
for five different halo mass intervals (as in i45i| ) . 

In order therefore the estimate the constants 60 and 
C2 within the context of the specific cosmological model 
used, we fit the DM halo bias results of the N-body simu- 
lation to the corresponding theoretical bias formula p.8p 
and find accurate fitting formulas for both parameters. 
These are 



(A.l) 



with = 0.857 ± 0.021 and /3 = 0.55 ± 0.06. and 



C2{Mh) = 



© 



(A.2) 



with = 1.105 ± 0.018 and ^ = 0.255 ± 0.005. In figure 
3 we show the simulation based values of the 69 and C2 
parameters as a function of DM halo mass together with 
the best fitted functions, provided in eq. (A.l) and (A.2). 

The cosmological dependence of these parameters is 
the subject of work in progress. 
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